Generalized Neighbor-Interaction Models Induced by Nonlinear Lattices. 
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It is shown that the tight-binding approximation of the nonhnear Schrodinger equation with a 
periodic hnear potential and periodic in space nonlinearity coefficient gives rise to a number of 
nonlinear lattices with complex, both linear and nonlinear, neighbor interactions. The obtained 
lattices present non-standard possibilities, among which we mention a quasi-linear regime, where 
the pulse dynamics obeys essentially the linear Schrodinger equation. We analyze the properties 
of such models both in connection with their modulational stability, as well as in regard to the 
existence and stability of their localized solitary wave solutions. 
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I. INTRODUCTION 

It is generally recognized that mapping of a nonlin- 
ear evolution problem, described by a partial differential 
equation, into a simplified lattice, representing a set of 
coupled ordinary differential equations, appears to be a 
useful tool either for numerical (or semi-analytical in the 
appropriate limits) study of the dynamics or for bring- 
ing intuitive understanding of the factors dominating the 
behavior of the systems. Examples of such approach 
are well known for a long time in solid state physics [!} 
(the description of an electron in a crystal in the tight- 
binding approximation), in optics ^] (the description of 
the electric field in arrays of waveguides), and more re- 
cently in the mean-field theory of Bosc-Einstcin conden- 
sates loaded in optical lattices (see e.g. [3, 4] for relevant 
reviews). In all mentioned cases, the periodicity is usu- 
ally associated with the linear properties of the system 
and the respective dynamics is approximately described 
by the discrete nonlinear Schrodinger (DNLS) equation. 

On the other hand, there has recently been an increas- 
ing interest in studying nonlinear models, where the non- 
linearity is also periodically modulated in space. Appli- 
cations of such models extend from the propagation of 
electromagnetic waves in stratified media Q to conden- 
sates of bosons and condensates of boson- fermion mix- 
tures in optical lattices. It turned out that the spa- 
tially dependent nonlinearity may dramatically change 
properties of the system, in particular the regions of exis- 
tence of coherent localized structures and especially their 
corresponding stability properties. 
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Natural questions that arise in this context are the 
mapping of the respective evolution equation into non- 
linear lattice and the description of the existence, sta- 
bility, and dynamics properties of such lattices. In the 
present paper we consider both of these issues. While 
the first issue is technical and can be straightforwardly 
addressed by means of the Wannier function expansion, 
as it was suggested in Q, the study of the properties 
of the emerging lattices is a much richer problem in the 
present setting; its richness stems from the fact that spa- 
tially dependent nonlinearity gives rise to complex non- 
linear inter-site interactions. Such additional forms of 
nonlinearity as the ones extracted below can significantly 
change the dynamical properties of the discrete system 
as it was shown in earlier research devoted to the spin 
waves in magneti c sy stems Q, electromagnetic waves in 
waveguide arrays [lOL and to applications to arrays of 
Bose-Einstein condensates [12] • In the above mentioned 
studies of nonlinear lattices, however, one common fea- 
ture was of crucial importance - that was the presence 
of dominant (or at least significant) on-site nonlinear- 
ity, which e.g., in the case of a BEG loaded in an op- 
tical lattice is typically about two orders of magnitude 
larger than the hopping nonlinearity. In the present pa- 
per, we systematically derive and consider a far more 
general class of lattice evolution equations, including the 
cases where the on-site nonlinearity is exactly zero. 

The organization of the paper is as follows. In Sec. |TT] 
we deduce the relevant lattice dynamical models with 
inter-site nonlinearity starting with the evolution equa- 
tion of the nonlinear Schrodinger (NLS) type with a spa- 
tially periodic potential and periodic nonlinearity. In Sec. 
mil we conduct the modulational stability analysis of the 
derived models. In sections llVl andlVl we examine the re- 
spective dynamical properties of the derived quasi-linear 
and nonlinear models, while Sec. IVII summarizes our find- 
ings and presents our conclusions, as well as directions of 
potential future interest. 
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II. MODEL EQUATIONS 
A. One-band approximation 

We start with the one-dimensional NLS equation 

i^ = -^+U{x)^ + gixM\^i;, (1) 

where and g{x) are the coordinate-dependent linear 
and nonlinear potentials respectively, both considered to 
be TT-periodic functions: U{x) ^ U{x + tt) and Q{x) = 
Q{x -I- tt) . We concentrate on the cases where the linear 
potential is an even function lA{x) — U{—x), while the 
nonlinearity may be either even or odd: G{x) = (tQ(—x) 
(hereafter cr = ±1). 

In order to map Eq. ([T]) into a lattice equation we 
follow [8] . To this end, we introduce the linear eigenvalue 
problem 

(]^ip^ {x\ 

^1^{x)Vaq{x) ^ £aq'^aq{x) (2) 

where ipaq{x) is a Bloch function, a > 1 and q stand 
for the band number and for the wavenumber in the first 
Brillouin zone: g G [—1, 1], and define the Wannier func- 
tions 



71 1/"^ 



(3) 



which constitute an orthonormal set of real and exponen- 
tially decaying functions ^1.3]. 

We seek the solution of Eq. ^ in the form of a series 



^j{x,t) = y^^Cna{t)Wna{x). 



(4) 



For the next consideration we notice that the Wannier 
functions of a-th band possess either even or odd parity 
WQa{x) = {—lY^°''WQa{—x) and are characterized by the 
property Wna{x) — woa{x — mr). 

Now we make the most crucial approximation of our 
model, namely that the continuum Eq. ([T]) can be accu- 
rately described within the one-band approximation. As 
it was shown in this assumption fails to describe the 
original continuous model when one studies dynamical 
processes associated with the generation of the frequen- 
cies belonging to the higher bands. However, it is reason- 
ably accurate in describing static solutions (in particular, 
localized modes) as well as their stability. Also, the lat- 
tices of generalized neighbor interactions derived below 
within the framework of the one-band approximation, are 
of interest in their own right, per their particularities and 
differences in comparison to other models of similar type; 

cf. Elm- 

We thus assume that only one band, say a-th one, is 
populated. Now, substituting the expansion (U) in Eq. 



(H)) we arrive at the equation (see [1, Q for more details) 

(5) 



"1,112, 13 



err wn"i"2n3 _ n 



where 



Q{x)Wna{x)Wmat{x)Wn2aAx)Wn-^aA^)dx (6) 



are the nonlinear overlap integrals, the coeffi- 

cients of the Fourier series expansion of the eigenvalue 
f ■ 

^aq • 

£aq = '^-e*""' . ^nc. = l f So.qe-'^^'^dq; (7) 

the overbar stands for complex conjugation, and an over- 
dot stands for the derivative with respect to time. In Eq. 
([5]) we have taken into account that in a general situation, 
for a periodic potential U{x) of rather large amplitude, 
the lowest bands are very narrow, and hence the Fourier 
coefficients Una decay rapidly with increasing n (so that 
I^^OqI ^ |wia| ^ |w2a|), which, in turn, allows us to ne- 
glect the coefficients tOna with n > 2. 

For the lowest bands the Wannier functions are well 
localized on the scale of one lattice period (and can be 
reasonably well approximated by the eigenstates of the 
linear oscillator) and thus the overlap integrals involving 
next-nearest neighbors (i.e., lattice minima separated by 
two lattice maxima) are negligibly small. This allows us 
to drop also the terms involving W^lalalat "^ith at least 
one pair of the upper indices satisfying \nj — nk\ > 2. 
We however emphasize, that it is of crucial importance 
to leave the nonlinear terms with hopping between the 
neighbor sites, which for specific choices of the nonlin- 
ear interactions Q{x) can be comparable with or even 
stronger than the on-site nonlinearity (see below). 

Now we use the symmetry of the integrals W with 
respect to permutations of the indices and introduce 



Wn = W"""" = w/oooo 



ytr TJ^n,n — l,n — l,n — 1 _TlAn,n,n,n— 1 



W9 = W",","~l,"-1 = w 



■W, 



1000 



1100 



(8) 



where a — 1 and a ^ —1 for G{x) even and odd and 
W,^ gix)wi^{x)wt-'{x)dx J = 0,1,2. (9) 



In the case of odd nonlinearity {a = —1) the terms Wq 
and W2 are always equal to zero due to the fact that the 
integrand in ([9]) is odd with respect to the points X = 
and X = 7r/2, correspondingly. We thus arrive at the 
equation 

iCn = UjQCn + LOliCn-l + Cn+l) + Wo\Cnf C„ 

+ crcn-icl + 2(7 1 

Cn I — 1 

-|-2|c„pC„+i +Cn+icl + Cr|c„+ipC„+i) 

+W2 (2|c„_i|2c„ + c„4_^ +c„4_^i + 2|c„+ipc„) ,(10) 
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where we have dropped the zone index a (e.g. c„^q. is 
redefined as c„, etc.) 

Eq. (fTO| is the main discrete model studied in the 
present paper. We notice that it has a Hamiltonian struc- 
ture: icn — dH / dcn with the Hamiltonian 



H 



E 



I l2 / - - N W^O, 

Wo|C„| + Wi (C„_1C„ + C„_iC„j + — |c„ 



+W^i(|c„_i|2 + (t| 

I ) {pn^n-- 



-2W^2|c„-inc„P 



^(c^c-Li+ctc^O 111) 



and with the standard Poisson brackets. Another inte- 
gral of motion is the sum N — ^„ |cnp, reflecting the 
conservation of the "number of atoms" (in keeping with 
the BEC motivation of our analysis) of the original Eq. 



Particular cases 



The deduced model (|10p allows for a number of inter- 
esting particular cases. First of all, however we notice 
that in the case of odd nonlinearity function, i.e., for 
(7 = — 1, neither purely even (c„ = c_„), nor purely odd 
(c„ — — c_„) solutions can exist. This follows directly 
from the symmetry of Eq. pO)) [or Eq. ([T])] . 

For illustration of the localized solutions we use the po- 
tential U{X) = — 3cos(2X) (both for the detailed calcu- 
lations of this section, and for the numerical simulations 
of the following sections). Also we restrict our consider- 
ations to the first band, i.e. we take a — 1. In that case, 
the linear overlap coefhcients for the first band are com- 
puted as Wo w -0.839, wi w -0.051, and uj2 ~ 0.002 

One can distinguish five cases as follows: 

Case 1: Wq = W\ = W2 = 0. This is a quasi-linear 
case, which is made possible, for example by an odd non- 
linearity (ct = — 1) of the form 



g{x) = sin(2a;) - 1.3706 sin(4a;) 



(12) 



(obviously a number of possible realizations of this and 
other cases reported below is naturally unlimited). Our 
approach both in this example and below is motivated by 
the nature of the lattice that we wish to construct (i.e., 
by the type of overlap integral that we wish to preserve 
or eliminate). For instance, in this example, the odd 
nonlinearity guarantees that Wq = W2 = 0, while the 
expression of Eq. IT^ uses one tunable parameter (the 
amplitude of the second harmonic) to achieve Wi = 0. 

Case 2: Wq = W2 = 0, Wi ^ 0, is achieved e.g. by 
choosing 

^;(a;) = 10sin(2a;) (13) 
((T = —1). Now the lattice model (ITU|) is reduced to 

iCn = t^0C,i + UJi{Cn-l + Cn+l) 



-Wi (|c„_ipc„_i - c„_ic^ - 2|c„pc„_i 



and in the case at hand Wi « 0.045. 

Alternatively, this case can be realized by choosing 
even nonlinearity ((7 = 1) 

Gix) = 0.0275 + 4.809 cos(2a;) - 10cos(4x). (15) 

Now Wi « 0.012 and the lattice equation reads: 

iCn = i^oCn + Wl(c,i-1 + Cn+l) 

+Wi (|c„_ipc„_i -f c„_ic^ -I- 2|c„|^c„_i 

+ 2|c„|^c„+i +c„+i4 + |c„+ipc„+i) . (16) 

Case 3: Wq ^ Wi ^ 0, W2 ^ is obtained, for in- 
stance, for the even nonlinearity (cr = 1) 

g{x) = -23.836 + 48. 882 cos(2j;) - 37.778 cos(4a;), (17) 



for which W2 ~ —0.0136. The model (|T0|) is now simpli- 
fied 

iCn = WoC„ -I- UJi{Cn-l + Cn+l) 

+W2 (2|c„_i|2c„ -Hc„4_^ -Hc„4^i + 2|c„+ipc„) .(18) 

Case 4- Wo — 0, 14^1,2 7^ can be achieved by using 
an even nonlinearity {a — 1) 



g{x) = 7.795 - 10 cos(2a;). 



(19) 



2|c„pc„+i 



i+lC„ - \Cn+l\ Cn+l) , 



(14) 



Now the overlap integrals are as follows Wi « 0.0148 and 
14^2 ~ 0.0045. The lattice model then reads 

iCn = LOQCn + U}l{Cn-l + Cn+l) 

+Wi (|c„_ipc„_i + c„_ic^ + 2|c„pc„_i 
-|-2|c„pc„+i + c„+ic^ + |c„+ipc„+i) 
+W2 (2|c„^i|2c„ + c„4_^ +c„4^i + 2|c„+ipc„) .(20) 

Case 5: |Wi|, IW2I < |Wo| is the standard case of the 
on-site nonlinearity (the DNLS equation): 

iCn = UjQCn + U!i{Cn-l + Cn+l) + Wo|c„pC„. (21) 

This form of the lattice dynamical model is obtained for 
generic nonlinearities, Ea. (|2T|) is well studied in the lit- 
erature (see e.g. [ij]) and that is why it will not be 
addressed in this paper. 



III. MODULATIONAL INSTABILITY 

As it is customary we start with the analysis of the 
modulational instability of plane-wave solutions of Eq. 

(for a general study of the modulational instability of 
the plane wave background in the DNLS-type equations 
see e.g. 

, [iE[i3)- This stability analysis is performed 
not only because it is of interest in its own right but 
also because the solitary wave solutions that we plan on 
constructing for the above presented lattices should be 
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produced on a dynamically stable background. Using 
the plane wave solution of the form: 



(22) 



where -F is a constant amplitude, we obtain the dispersion 
relation (in the absence of the previously considered cubic 
onsite terms) 

+ 2uji cos(g) + iWiF'^ia + 1) cos(g) + 
2W2F^[2 + cos{2q)]. (23) 

To examine the linear stability, we perturb the plane 
wave solution in the form: 

Cn = {F + ^e*(<^"-"*) + ije-'^*?"-"*) )e'(«"-"*) , (24) 

with |A| , <C \F\ and linearize with respect to A and 
B. As a result we obtain two branches of the linear ex- 
citations fll.2{Q)- 



(25) 



where 

D = 

M_ = 
M+ = 
A = 



ijj -ijjQ- AWi{l + o-)F^ cos(g) - 

AW2F^ (cos(Q) + 1) , 

-2 {uji + 2Wi{l + a)F'^) sin(Q) sin(g) - 

414^2^2 sin(Q)sin(2g), 

2 (wi + 214^1 (1 + a)F^) cos(Q) cos(g) + 

AW2F'^ cos(g) cos(2g) , 

A[Wi{l + (t){1 + cos(g)) cos(q)+ 

M^2(2cos(g) +cos(2g))]V 

4(1 - af svn^{q)Wl{l - cos{Q)f. 



Let us consider different special cases for cr — —1 (the 
case (T = 1 is investigated in [loj). 

a. Homogeneous background is described hy q — Q. 
Now the two branches of the solutions collapse and the 
dispersion relation acquires the form 



16(^1 + 214^2^2) sin^ 



(Wi 



61^2^2) sin^ 1^1" 'iW2F' 



(26) 



The carrier field is stable if and only if the two following 
conditions are satisfied 



{loi + 2W2F^) W2 < 0, 
{uji + 2W2F^) {uji + 3W2F^) > 0, 



(27) 



where the first of these conditions demands non- 
negativity of the coefficient k in the expansion of Eq. ([26]) 
of the form of fl^{Q) — kQ^, which is valid for small Q. 
The second condition demands non- negativity of fl^{Q) 
at the zone boundary, Q = tt. 



We point out that the long wavelength excitations' 
group velocity dispersion is given by 



— ^2W2[2W2F^ + {8W2-6 + LUi)F^+Auji] 

Q=0 

and takes zero values for 



(28) 
1/2 



4W^ ^ 



'6-LU1-8W2 
. 4W^ 



2u;i 



This last condition determines the domain of the param- 
eters where shock waves can be observed [18{ . 

b. Staggered background corresponds to g = tt. The 
dispersion relation is 



16(wi - 2W2F^)sm-' ^ 



(uji - 6W2F^) sin^ ( ^ ) + 31^2^^ 



(29) 



Similar to (|27p we introduce stability criteria in the form 

(wi - 2W2F^) W2 > 0, 
(wi - 2W2F2) (wi - 3W2F2) > 0. (30) 

Note that the dispersion relations for q = 0, tt do not 
depend on the coefficient Wi. 

c. Phase alternating background where q — tt/2. 
The dispersion relation is 

fli.2 ~ —2uji sin(Q) ± 
4V2F^ sin (^^^ ^JiWi + W^)cos{Q)-Wl (31) 

The condition for the presence of modulational instability 
in this case is 



W? 



(32) 



IV. QUASI-LINEAR MODEL (CASE 1) 

Let us now turn to the quasi- linear model of Case 1. 
We notice that the term "quasi-linear" is used in order 
to emphasize the existence of solely higher order contri- 
butions from more distant neighbors which are small but 
not necessarily zero. This provides us with an excellent 
benchmark of our derivation since in the discrete linear 
case, the dynamical equation can be solved explicitly and 
subsequently compared to the full results of the original 
partial differential equation (for which the discrete model 
was developed as an approximation). In particular, it 
is well-known that the linear discrete case, with a com- 
pactum of initial data c„ (0) — ASno has a solution of the 
form 



Cn{t) = A (-i)" exp(-iwot) J„(2wit) 



(33) 
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where J„ is the Bessel function of order n. 

We have tested this analytical prediction of the dis- 
crete model in the partial differential equation ([1]) with 
the "nonlinear potential" of Eq. . The results of our 
numerical simulations can be found in Fig. [H which high- 
lights the excellent agreement between the analytical and 
the numerical results. This also serves to showcase the 
accuracy of the reduction via the tight-binding approx- 
imation of the original partial differential equation. We 
will hereafter focus more on the details of these discrete 
models and of their solitary wave solutions. 

The remarkable accuracy of the tight-binding model 
in the case at hand can be easily understood. Indeed, 
we are dealing with a perfect lattice (i.e. having no de- 
fects). Inter-band transitions, which are the cause of the 
failure of the one-band approximation when they exist, 
are only due to the nonlinear coupling of bands, and 
exactly this factor is anomalously small for the chosen 
nonlinearity. The first indication on this fact is given 
by the zero contributions of Woa,2- Next, due to the 
symmetry of the Wannier functions, and the symmetry 
of the nonlinearity, one concludes that W^^Jp = for 
all a and /3, i.e., there exists no tunneling between the 
same sites of two different bands. Even more generally 
(also due to the symmetry), W^"'^^™'+2^™'" = (for ar- 
bitrary integers /3, 7, n and m). For completeness we 
have checked the numerical values of the other inter-band 
overlap integrals for the three lowest bands. The integrals 
greater than 0.01 are as follows (the ones obtained by the 
symmetry reductions are not shown) VF^m — 0.0465, 
W^222i° = -0.1360, H/2m° = -0.0228, W^^^^ = -0.0432, 
= -0.0301, Wl^^^ = -0.0117. The coefHcients 
describing energy transfer from the first to the second 
and the third bands are W^f^f = 0.0465, W^^f^ = 0, 
W^2m = 0.0077, Wl^^^ = -0.0432, W^l^f = 0.0004, 
VFgYn — 0.0011, i.e., either have relatively small value or 
are identically zero. This explains the high accuracy of 
the one-band approximation. 

Finally we notice, that in the case at hand the disper- 
sion relation (PS)) takes the form independent of the wave 
amplitude F 

n = ±2wi [cos (q) - cos (Q T q)] , (34) 

which is the dispersion relation for linear phonons, and 
hence no instabilities can occur. 



V. NONLINEAR MODELS (CASES 2 AND 3) 

The above confirmation of the high accuracy of the 
tight-binding model for the description of the pulse dy- 
namics motivated us to study the coherent structure so- 
lutions of Eq. (fTU)) and their properties for the Cases 2 
and 3 singled out in Sec. IIIBI It is convenient to join the 
equations of the Cases 2 and 3 and to consider for a = 1 



c 



20 40 60 80 100 
t 




t 



FIG. 1: The top and middle panel show the space-time con- 
tour plot of the solution in the PDE of Eq. ^ [top panel] and 
the discrete model analytical prediction (|33|) [middle panel]. 
In the latter case, space is normalized over the period of the 
linear potential, so that it can be compared with the lattice 
results. To accentuate the excellent agreement between ana- 
lytical and numerical results, the bottom panel shows the time 
evolution of the amplitude at the central site of the configura- 
tion compared between the PDE numerical result (solid line) 
and the discrete equation analytical result of (|33} (dashed 
line). The two are practically indistinguishable. 
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W=0 (Case 2) 
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FIG. 2: Equilibrium soliton solutions found in the Case 2 
{W2 = 0): (a) pulse, (b) kink, and (c) anti-dark structure with 
oscillatory wave structure around the background, hereafter 
termed wave. On-site structures are presented here, but inter- 
site ones can also be constructed. Solutions in (a) and (b) 
can be stable but we were unable to find a stable structure of 
the form presented in (c) (see Sec. IVC[l . Model parameters 
corresponding to panels (a) to (c) are depicted by dots in 
Fig. [3] marked by the capital letters A to C, respectively. 
Parameters are: (a) Wi = -0.0112, lo = -0.99, (b) Wi = 
0.0112, u = -0.9, (c) Wi = 0.0112, uj = -0.75. 

the model 

ic„ = IVoCn + wi(c„_i -t- C„+i) 
+ Wi {\Cn^i\^Cn-l+Cn-lcl+2\Cnfcn-l 
+ 2\Cnfcn+l +Cn+lCn + |Cn+l|^C„+i) 
+W2 (2|c„_i|^c„ + c„4_^ + C„C^+i + 2\Cn+lfcn) -(35) 

Let us seek stationary solutions of Eq. ((35| of the form 

cn (t) = fne-"^\ (36) 

with real /„. Using this ansatz, we obtain the equation 
for the amplitudes 

(UJO - Uj) f„ + UJl ifn-l + f„+l) (37) 

+W, + -I- 3/2 (/„+i + /„_i)] 

+ 3W2fn {fLl + fn+l) = 0. 

We attempt to find the pulse (bright soliton) and 
kink (dark soliton) solutions. Our strategy in search- 
ing for these solutions will be as follows. We will first 
formulate the necessary conditions for the existence of 
the soliton solutions considering their carrier constant- 
amplitude field and also the soliton tail solutions. This 



FIG. 3: Parameter plane (Wi,uj) in the Case 2 {W2 = 0) with 
indicated regions of existence of three types of solitons shown 
in Fig. [2] and their staggered analogues. The pulse solution 
exists for Wi < and lj < u>t, while the staggered pulse for 
Wi < and w > Ua- The kink and staggered wave can exist 
for Wi > and < < Uc, while the staggered kink and 
the wave for Wi > and Ud < ui < Ua- The values of special 
frequencies Ua to ujd are given by Eq. (|46p and Eq. (|47p . 



will narrow the domain of parameters where such so- 
lutions can be expected. Then, we will attempt to 
construct the desired soliton solutions and subsequently 
study their stability. 

The staggered and non-staggered stationary solutions 
are connected by the following symmetry relation p^ : 
if /„ is a solution of Eq. ([57]) for definite Wi, W2, and 
UJ < ujo {uj > ujq), then (—!)"/„ is a solution for Wi = 
Wi, W2 = —W2, and UJ — 2ujq — uj > ujq {uj < ujq). 
We also note that the stability analysis of any stationary 
solution can also be done, without loss of generality, for 
only, say, non-staggered carrier field. This is so because 
the dynamics in the vicinity of the stationary solution is 
governed by Eq. (|5ip (see below) which is invariant with 
respect to the following transformation: e„ — > (— l)"e„. 
In ^ (-l)"/n, uj^2ujQ- UJ, W2 ^ ~W2, and t ^ -t. 
Bearing this in mind, in the following we will discuss only 
stationary solutions with a non-staggered background. 



A. Constant amplitude solution and soliton tails 

A pulse solution must satisfy the boundary conditions 
|cnp — > for n ±00, while for the kink solution one 
must have |c„p ^ > for n ±00. Thus, the 
existence and stability of the carrying field solution, de- 
scribed by formula with real F, is a necessary con- 
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FIG. 4: Equilibrium soliton solutions obtained in Case 3 with 
Wi = 0: (a) pulse, (b) kink, and (c) kink with oscillatory 
tail, called oscillatory kink. On-site structures are presented 
here, but inter-site ones can also be constructed. All these 
solutions can be stable, as it will be shown in Sec. IV CI Model 
parameters corresponding to panels (a) to (c) are depicted 
by dots in Fig. [5] marked by the capital Latin letters A to 
C, respectively. Parameters are: (a) W2 ~ —0.0136, ui = 
-0.9704, (b) W2 = 0.0136, uj = -0.9, (c) W2 = 0.0136, u = 
-0.85. 



FIG. 5: Parameter plane {W2,lo) in Case 3 {Wi — 0) with 
indicated regions of existence of three types of solitons shown 
in Fig. |4]and their staggered analogues. The pulse exists for 
W2 < and UJ < uji,, while the staggered pulse for W2 > and 
u) > u!a- The kink exists for W2 > and ujb < < ujf, while 
the staggered kink for W2 < and uie < uj < u)a- Finally, the 
oscillatory kink exists for W2 > and ujf < uj < uja, while 
the oscillatory staggered kink for W2 < and luq < uj < uje. 
Values of special frequencies are given by Eq. (|46p and Eq. 

(BHl). 



dition for the existence and stability of the sohton sohi- 
tions. 

One always has the trivial solution F ^ and from 
the expression (|23p one can have two nonzero solutions 
with 



(38) 



8W1+6W2 ' 

if the expression in the right-hand side of Eq. is 
positive. In the case 8Wi + 6W2 = 0, F can be arbitrary 
if w = + 2ti;i, but we will not study this very special 
case. 



Substituting Eq. ([381) into Eq. we obtain the 

spectrum of the carrier field with F^ > 0, whose stability 
criteria (in full analogy with Sec. IIII[) are 

— LOo + 2aJi) + W2{lo — ujq + uJi)] x 
(41^1 -I- 3W2) {lu-luo~ 2lui) < 0, 

[Wl{uJ - UJO + 2uJl) + W2{iL! - OJQ + OJl)] X 

[2Wi{lj -LUO + 2wi) + 31^2 - Luo)] > 0. (39) 

From the asymptotic properties mentioned above, one 
can express the soliton tails as 



(40) 



at \n\ 00, where small ^„ are real and are independent 
on t. Substituting Eq. (|40)) into Eq. (j37p and linearizing 
with respect to ^„ one obtains 



with 



(3^LUo-LU + l2WiF^ + QW2F^, 
7 = wi -I- 61^1^2 + 61^2 i^^. 



(41) 



(42) 



Seeking solutions to Eq. (|41]l in the form 

^„ ~ C±z", n ±00 
with complex z, we come to the characteristic equation 
7^2 + /32 + = 0. (43) 
Thus, z is one of the roots zi^2 



1 



zi = 




27 y 47 



(44) 



providing |^„| — > as |n| — > 00. 

For the soliton, which is either a homoclinic or a hete- 
roclinic of the map, generated by Eq. ((57)) . ^„ = must 
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FIG. 6: Pulse in Case 2 (W2 — 0). (a) Space-time evolu- 
tion of |c„(t)p showing the dynamics of the pulse placed at 
t = Q asymmetrically with respect to the lattice. The pulse 
undergoes periodic oscillations in the vicinity of the stable 
inter-site configuration, (b) Time variation of the norm of the 
four central particles. Parameters: Wi — —0.012, W2 = 0, 
uj — —0.99, which corresponds to the point A in Fig. [S] 

be a hyperbolic point which happens only if the roots zi^2 
are real, i.e., if 

^32 _ 4^2 ^ Q (-45^ 

The sign of z specifies the type of the tail solution: 
z > corresponds to the tail decaying monotonically 
with distance from the soliton's center, while z < means 
that the decaying tail solution oscillates near the carrier 
solution F. 

The absolute value of z characterizes the degree of lo- 
calization of the tail. If [zj is small, then the tail solution 
is weakly localized, otherwise it is strongly localized. We 
found that, in many cases, there is a correlation between 
the degree of localization of the tail solution and that of 
the soliton itself. 

The tail solution was found from the linearized equa- 
tion and it can only provide necessary conditions for the 
existence of a soliton solution. The nonlinear terms, de- 
pending on their structure and the values of the corre- 
sponding coefficients, can either lead to unbounded so- 
lutions or to non-localized solutions, or to the desired 
bounded and localized soliton solutions. 

Having a tail solution one can attempt to construct the 
corresponding soliton solution. To do so, we use the tail 
solution defined by Eq. ([^0]) to set the initial values for 
fn-i and /„ and then find fn+i fi'om Eq. (|37p . solving 
the cubic (for W2 — 0) or the quadratic (for Wi — 0) 
algebraic equation. The proper choice of the integration 
constants C± systematically allows one to obtain equilib- 
rium on-site or inter-site soliton solutions, if they exist 




FIG. 7: Pulse in Case 3 {Wi — 0). (a) Space-time evolution 
of |c,i(f)p showing the dynamics of the pulse placed at t = 
asymmetrically with respect to the lattice. The pulse under- 
goes periodic oscillations in the vicinity of the stable on-site 
configuration, (b) Time variation of the norm of the central 
particle and its two nearest neighbors. Parameters: Wi — 0, 
W2 = —0.0136, Lu — —0.9704, which corresponds to the point 
A in Fig. [3 

To conclude this section we summarize the necessary 
conditions for the existence of pulses and kinks. 

Pulse in Cases 2 and 3. Since the carrying field for 
pulses with F — Q always exists, there remains only one 
necessary condition, namely, the condition of the exis- 
tence of the tail solution of Eq. (|45p. For F = this 
condition is satisfied for any Wi and W2 and for both 
Cases 2 and 3 when 

w > Wo - 2wi = -0.73732741 = cja, 

w < Wo + 2wi = -0.94041721 = uj^- ^ ' 

These conditions simply state that the localized pulses 
must be located outside the phonon band of the spectrum 
given by the interval [wb, Wq]. 

The necessary conditions of the existence of a kink in- 
clude the condition of the existence of the carrier field 
with > 0, Eq. the stabihty condition for the 

carrier field, Eq. (|39p . and condition Eq. (|45p of the 
existence of the tail solution. 

Kink in Case 2. All three necessary conditions are 
satisfied when Wi > and LUh < lu < cua, and they are 
not satisfied simultaneously for Wi < 0. 

Kink in Case 3. All three necessary conditions are 
satisfied when W2 > and ojb < uj < ujq, while for W2 < 
they are satisfied for loq < uj < uia- 

We also specify the frequencies at which z changes sign. 
In the Case 2 the frequencies are 

ujo - (2/3)wi = -0.80502401 = ujc, 

Wo -t- (2/3 Vi = -0.87272061 = Wd, ^ ' 
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FIG. 8: Kink in Case 2 (W2 = 0). (a) Space-time evolu- 
tion of |c„(t)p showing the dynamics of the kink placed at 
t — asymmetrically with respect to the lattice. The kink 
undergoes periodic oscillations in the vicinity of the stable 
on-site configuration, (b) Time variation of the norm of the 
central particle and its two nearest neighbors. Parameters: 
Wi = 0.012, W2 ^ 0, to = -0.9, which corresponds to the 
point B in Fig. [S] 

while in the Case 3 they are 

t^o-tc^i = -0.78809986 = We, /-J.^ 
Wo + t^i = -0.88964476 = w/. 

B. Soliton solutions 

case of It is well-known [l^ [2l[ that the standard 
DNLS equation (|2ip possesses multiple branches of local- 
ized solutions which are parametrized by the frequency 
detuning outwards the phonon band. Multiple branches 
of the stationary localized solutions were also obtained 
in a different model in the presence of both linear and 
nonlinear latices Q. This allows one to conjecture that 
any of the lattices introduced in Sec. should also pos- 
sess multiple branches of the localized solutions and kinks 
(discrete dark solitons). Although a thorough study of 
each of the cases is by itself a nontrivial problem that 
requires lengthy considerations, to present a panoramic 
view of the possible nonlinear modes in the above lat- 
tices, in the present paper we restrict our study to some 
representative examples, referring, in most cases, to the 
lowest branches. 

Case 2 (W2 = 0). 

For the nonlinearity of the Case 2, using the method 
described in Sec. IV Al we obtain examples of local- 
ized solutions of the Eq. namely: pulse, staggered 
pulse, kink, staggered kink, wave (an anti-dark struc- 
ture in the form of an oscillatory wave on a non-zero 



FIG. 9: Kink with oscillatory tail in Case 3 {Wi = 0). (a) 
Space-time evolution of |c„(f)p showing the dynamics of the 
kink with oscillatory tail placed at t = asymmetrically with 
respect to the lattice. The kink with oscillatory tail under- 
goes periodic oscillations in the vicinity of the stable on-site 
configuration, (b) Time variation of the norm of the central 
particle and its two nearest neighbors. Parameters: Wi — 0, 
W2 = 0.0136, Lj = -0.88. The kink with oscillatory tail has 
frequency close to ujf — —0.8896 and it becomes unstable far 
from this line (see Fig. llSp . 



background), and staggered wave. The non-staggered 
solutions are presented in Fig. [2] for such a parame- 
ters: (a) Wi = -0.0112, Lu = -0.99, (b) Wi = 0.0112, 
Lu = -0.9, (c) Wi = 0.0112, u = -0.75. The corre- 
sponding staggered solutions can be constructed using 
the staggering transformation, described above. On-site 
equilibrium structures are shown but one can also obtain 
the inter-site ones. 

In Fig. O the regions of existence of various solutions 
are indicated on the parameter plane [Wi^oj). Solutions 
presented in Fig. [5] (a) to (c) have parameters shown 
by dots marked by the capital letters A to C, respec- 
tively. Recall that pulses can exist in the two frequency 
ranges, uj < ujjj and w > uja, for any Wi. However, they 
were found only for Wi < 0, while for Wi > the itera- 
tions initiated by the tail solution resulted in unbounded 
structures. In the portion of the plane with Wi > and 
Wf, < a; < Wc kinks [see Fig. [5] (b)] and staggered waves 
were obtained. On the other hand, in the portion with 
Wi > and tOd < tu < LUa, we could construct staggered 
kinks and waves [see panel (c) of Fig. [2]. It is important 
to note that j^j is close to 1 near the lines u = uja and 
uj = ujb, where pulses, kinks, and waves were found to 
be wide; z diverges (or vanishes) at oj — ujc for staggered 
carrier field and at lu — Ud for non-staggered carrier field; 
in the range of uJd < uj < ujc, z is always negative and 
close to 0, resulting in rapidly oscillating or sharply lo- 
calized (and typically unstable) solutions found from the 
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tail construction. 
Case 3 {Wi = 0). 

In the Case 3, three examples of the non-staggered 
stationary soliton solutions presented in Fig. |4] were 
found. Shown are: (a) pulse, (b) kink, and (c) kink 
with oscillatory tail, referred to as oscillatory kink. On- 
site structures are presented here, but inter-site ones can 
also be constructed. Model parameters corresponding 
to panels (a) to (c) of Fig. |4] are depicted by dots in 
Fig. [5] marked by the capital letters A to C, respec- 
tively. Parameters for the non-staggered solutions are: 
(a) W2 = -0.0136, uj = -0.9704, (b) W2 = 0.0136, 
LO = -0.9, (c) W2 = 0.0136, UJ = -0.85, while the cor- 
responding staggered solutions, as in the Case 2, can be 
obtained using the staggering transformation. 

The pulse tail solution (with 2; > 0) exists for lo < ujh 
but the pulse itself exists in this region only for W2 < 0, 
while for positive W2 we obtained unbounded solutions. 
Similarly, the staggered pulse tail solution (with z < 0) 
exists for uj > uja but the staggered pulse itself exists in 
this region only for W2 > 0, while negative W2 leads to 
unbounded solutions. 

The kink exists for W2 > and lui, < uj < ujf, while the 
staggered kink for W2 < and uJe < cu < oja- Finally, the 
oscillatory kink exists for W2 > and lo f < uj < luq, while 
oscillatory staggered kink for W2 < and ujq < uj < uje- 

In Fig. [5] and Fig. [7] we contrast the behavior of pules 
in Cases 2 and 3, respectively. In both figures (a) shows 
the space-time evolution of |c„(i)p, while (b) shows time 
variation of the norm of the central particles. On pur- 
pose, we did not optimize the choice of C± to get equi- 
librium on-site or inter-site solutions. As a result, the 
ensuing profiles are non-stationary and, due to the pres- 
ence of the Peierls-Nabarro potential, they oscillate in 
the vicinity of stable configurations. It turns out that in 
the Case 2 (Case 3) the inter-site (on-site) configuration 
is stable. This conclusion will be confirmed in Sec. IV CI 
This illustrates the interes ting phenomenon of potential 
inversion of stability (cf. [l3|) in comparison with the 
standard DNLS mode Parameters in Fig. [5] are: 

Wi = -0.012, W2 = 0,uj = -0.99, which corresponds to 
the point A in Fig. [31 Parameters in Fig. [7] are: Wi — 0, 
W2 = —0.0136, UJ — —0.9704, which corresponds to the 
point A in Fig. O 

Similar results for dark solitons are presented in Fig. 
[5] (kink in Case 2) and Fig. [5] (kink with oscillatory tail 
in Case 3). One can see that in both cases the on-site 
structures are stable and this will be confirmed in Sec. 
rVCl Parameters in Fig. M are: Wi = 0.012, W2 = 0, 
UJ = —0.9, which corresponds to the point C in Fig. [31 
Parameters in Fig. M are: Wi = 0, 1^2 = 0.0136, uj = 
-0.88. 
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FIG. 10: Spectra of (a) unstable on-site and (b) stable inter- 
site pulses. The parameters (Case 2): Wi = —0.012, W2 — 0, 
UJ = —0.99 correspond to the point A in Fig. [3] 
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FIG. 11: Spectra of (a) stable on-site and (b) unstable inter- 
site kinks. The parameters (Case 2): Wi ^ 0.012, W2 = 0, 
u! — —0.9 correspond to the point B in Fig. [3] 



IV Bl We consider the following perturbed form of the 
solutions. 



^{t) = [/„ + e„(i)]e- 



(49) 
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C. Stability of soliton solutions 

Let us now study the stability of the stationary soli- 
ton solutions of the form of Eq. (|351) described in Sec. 



FIG. 12: Spectra of (a) unstable on-site and (b) unstable 
inter-site wave. The parameters (Case 2): Wi — 0.012, W2 = 
0, a; = —0.75 correspond to the point C in Fig. [31 
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where the small complex perturbation e„(t) is expressed 
as follows, 



(50) 



Substituting Eq. ^ into Eq. ^ we find that e„(i) is 
governed by the following linearized equation, 



ie„ = (wo -t^)e„ +tJi(e„_i +e„+i) 



+Wi 



+W2 



/n-l^n-l + 2/„(/n-l + fn+ljZn 



2/n/n-l^n-l + (/n-1 + /n+l)-^n 
+ 2fnfn+lZn+l , 



(51) 



where Z„ = 3a„ + «6n- Then, separating real and imagi- 
nary parts of Eq. (j5ip we derive the following system. 



K \ / b 
J j I a 



(52) 



where vectors a and b contain a„ and 6„, respectively, 
while the nonzero coefficients of matrices K and J are 
given by, 

K„,„_i = -aji-3Wi(/^i+/2)-6W^2/„-i/„, 

Kn.n = W - t^O - 6VKi/„(/„_i + /„+i) 

- 3M^2(/^i+/^+i), 

i^„,„+l = -CU,-3Wi{f^+f^+,)~6W2fnfu+l, 

J„,„_l = Cc^i + M^l(/^i + /,^) + 2W^2/„-l/n, 

Jn,™ = ^0 - + 2W^l/„(/„-l + /n + l) 

+ W^2(/^l+/^+l), 



n+1- 



(53) 



In the above expressions, n = l,...,Af, where Af is the 
number of lattice points. For pulses and kinks, we used 
periodic and anti-periodic boundary conditions, respec- 
tively. 

A stationary solution is characterized as linearly stable 
if and only if the eigenvalue problem 



K \ / b 
J j I a 
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(54) 



results in nonpositive real parts of all eigenvalues 7. 

The results of the stability analysis for the equilibrium 
structures reported in Sec. IV Bl are presented in Figs. 
[To] to [12] for the Case 2 and in Figs. [13] to [15] for the 
Case 3. The presented spectra contain (i) the vibration 
frequencies of the homogeneous background given by Eq. 
P5| ; (ii) a pair of zero- frequency modes corresponding to 
the invariance with respect to the phase shift; (iii) they 
also can include soliton internal modes falling outside 
the phonon band, see, e.g.. Fig. [TJ] (a). As was already 
mentioned, the spectra of stable structures do not possess 



eigenvalues with positive real parts, while those of the 
unstable ones have such eigenvalues. Now we turn to the 
discussion and comparison of the spectra in the Cases 2 
and 3. 

Case 2 (W2 ^ 0). 

Spectra of the on-site and inter-site pulses are pre- 
sented in Fig. [Tn](a) and (b), respectively. Interestingly, 
the inter-site configuration is stable while the on-site one 
is unstable. This type of instability is typical for the dis- 
crete system with Peierls-Nabarro potential, although in 
the standard cubic onsite nonlinearity case, the results 
are entirely reversed in comparison to the present case 
[1^ (e.g., the on-site pulse is stable, while the inter-site 
features a real eigenvalue pair). This indicates that the 
Case 2 nonlinearity results in the shape of the Peierls- 
Nabarro potential having a maximum (minimum) for the 
on-site (inter-site) pulses. Notice that as discussed in 
[l0| . such inversions of stability may occur in such gen- 
eralized models, upon varying their relevant parameters 
(such as Wi in the present case). We will see that for 
the Case 3 nonlinearity the situation for the pulse is ex- 
actly the opposite. The parameters used in this case are 
Wi = -0.012, W2 = 0,uj = -0.97, corresponding to the 
point A in Fig. [3] 

FigurefTTjshows results for the kink structures: the on- 
site kink in (a) is stable while the inter-site one in (b) is 
unstable. Here the location of maxima and minima of the 
Peierls-Nabarro potential is the same as in the classical 
discretization. The parameters Wi = 0.012, W2 = 0.0, 
uj = —0.9 correspond to the point B in Fig. [3] 

Finally, in Fig. [12] we show that (a) the on-site wave 
and (b) the inter-site wave are both unstable. The pa- 
rameters Wi — 0.012, W2 = 0, and w = —0.75 corre- 
spond to the point C in Fig. [3] 

Case 3 (Wi = 0). 

Results of the stability analysis are presented in Figs. 
[T^KTSl for the three soliton solutions displayed in panels 
(a) to (c) of Fig. [U respectively. The left panels show the 
spectra of the on-site structures and the right panels show 
the same for the corresponding inter-site structures. One 
can see that, in contrast to the Case 2, where the inter- 
site pulse was found to be stable, in the Case 3 the inter- 
site structures are always unstable (this is analogous to 
the case of the standard cubic discrete model with the 
on-site nonlinearity). This indicates that in the Case 3, 
the on-site (inter-site) structures are situated in the wells 
(on the peaks) of the Peierls-Nabarro potential. On the 
other hand, panels (a) in Fig. [T3] to Fig. [T^] present 
purely imaginary spectra for the on-site configurations, 
and this indicates that, for the corresponding values of 
model parameters, all three types of equilibrium solutions 
are stable. However, in Fig. [15] we demonstrate that the 
on-site configuration of the oscillatory kink is stable at 
uj = 0.88 [see panel (a)] but it can become unstable e.g. 
at oj = 0.85 [see panel (c)] with other parameters being 
unchanged and equal to Wi = 0, W2 = 0.0136. 
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FIG. 13: Spectra of (a) stable on-site and (b) unstable inter- 
site pulses. Oscillatory motion of the pulse in the vicinity of 
the on-site configuration is shown in Fig. [G] The parameters 
(Case 3) Wi ^ 0, W2 = -0.0136, uj = -0.9704 correspond to 
the point A in Fig. O 
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FIG. 15: The top panels show spectra of (a) stable on-site 
and (b) unstable inter-site oscillatory kinks with uj = —0.88. 
The bottom panels show the same for lu = —0.85 when both 
on-site and inter-site oscillatory kinks become unstable. Note 
that the type of instability in (c) is different from that in (d) . 
In (c) there are two pairs of complex-conjugate eigenvalues, 
while in (d) there is a pair of purely real eigenvalues. The rest 
of the parameters are (Case 3) Wi = 0, W2 = 0.0136. 



FIG. 14: Spectra of (a) stable on-site and (b) unstable inter- 
site kinks. The parameter values (Case 3) Wi — 0, W2 = 
0.0136, u — —0.9 correspond to the point B in Fig. [5] 



VI. CONCLUSIONS 

In this work, we have illustrated the potential for for- 
mulation of a rich variety of tight-binding nonlinear lat- 
tice dynamical models, stemming from the complex in- 
terplay of linear and nonlinear periodic potentials in the 
Gross-Pitaevskii equation. We have examined various 
particular possibilities, including quasi-linear models, as 
well as strongly nonlinear models where the nature of the 
coupling between the neighbors is itself nonlinear. Fur- 
thermore, in the nonlinear ones among our models we 
have studied the potential for the existence and stabil- 
ity, as well as the dynamics of localized solutions. More 
specifically, we have reported that the discrete model of 
Eq. dSni) with Wi^0,W2=0 (Case 2), and Wi = and 
W2 (Case 3) supports a number of localized station- 
ary solutions depicted in Fig. Oand Fig. S) Interestingly, 
the nonlinearity of the Case 2 results in stable inter-site 
pulse and stable on-site kink (both in staggered and non- 
staggered forms). On the other hand, the nonlinearity 
of the Case 3 supports only on-site stable localized so- 
lutions of three different types, namely, pulses, kinks. 



and kinks with oscillatory tails (all in both staggered and 
non-staggered forms). These results suggest that in the 
discrete models with nonlinear terms including interac- 
tions between nearest neighbors the profile of the Peierls- 
Nabarro potential can change qualitatively depending on 
the structure of the nonlinear terms and on the type of 
the coherent structure. In fact, it has been demonstrated 
that there exists a wide class of non-integrable discrete 
models of this sort where the Peierls-Nabarro potential is 
precisely equal to zero and equilibrium stationary solu- 
tions can be placed anywhere with respect to the lattice 
points [2^. 

A very natural extension of the present work would be 
to consider similar types of reductions in higher dimen- 
sional settings and to examine in detail the particular 
localized solutions that may emerge in the resulting lat- 
tice models. In particular, higher dimensionality may of- 
fer the potential for solutions with topological charge; it 
would therefore be relevant to examine under what con- 
ditions such solutions may be stable and how the relevant 
results relate to the original continuum model. 
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